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The coupled equations for the scalar modes of the linearized Einstein equations around 
Schwarzschild 's spacetime were reduced by Zerilli to a 1+1 wave equation d 2 ^/ z /dt 2 + TL*S> z = 0, 
where TL = —d 2 /dx 2 + V(x) is the Zerilli "Hamiltonian" , x the tortoise radial coordinate. From 
its definition, for smooth metric perturbations the field ^ z is singular at r s — —6M/(£ — 1)(£ + 2), 
with £ the mode harmonic number. The equation ty z obeys is also singular, since V has a second 
order pole at r 3 . This is irrelevant to the black hole exterior stability problem, where r > 2M > 0, 
and r s < 0, but it introduces a non trivial problem in the naked singular case where M < 0, then 
r 3 > 0, and the singularity appears in the relevant range of r ( < r < oo). We solve this problem 
by developing a new approach to the evolution of the even mode, based on a new gauge invariant 
function, that is a regular function of the metric perturbation for any value of M. The relation 
O | of to ^ z is provided by an intertwiner operator. The spatial pieces of the 1 + 1 wave equations 

D ■ that and \& z obey are related as a supersymmetric pair of quantum hamiltonians TL and TL. For 

M < 0, TL has a regular potential and a unique self-adjoint extension in a domain T> defined by a 
physically motivated boundary condition at r = 0. This allows to address the issue of evolution of 
gravitational perturbations in this non globally hyperbolic background. This formulation is used to 
complete the proof of the linear instability of the Schwarzschild naked singularity, by showing that 
a previously found unstable mode belongs to a complete basis of TL in D, and thus is excitable by 
generic initial data. This is further illustrated by numerically solving the linearized equations for 
suitably chosen initial data. 
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I. INTRODUCTION 



> 
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The linear stability under gravitational perturbations of the negative mass Schwarzschild spacetime was first consid- 
ered in where a proof of stability for the vector (or odd) modes is given. For the scalar (even) modes, reconsidered 
in 0, the problem is far more subtle, because the behaviour of the Zerilli potential V z [H, 0] at x — (which 
0^ ' corresponds to the r — Schwarzschild singularity) implies a one parameter ambiguity [l[ in boundary conditions 
at this point (parameterized by 9 € S 1 , see equation (j2"T) below), and also because V z has a second order pole at 
r = r s := — 6M /(£ — 1)(£ + 2) which falls within the domain of interest for negative M. None of these problems are 
present in the positive mass case, for which the relevant range is r > 2M (mapped over — oo < x < oo), and r s < 0. 

The ambiguity in boundary conditions at x = was addressed to in [3,12], where it was shown that ^ z ~ x 1 ! 2 as 
x — > + is to be selected in order that the first order corrections to the Riemann tensor algebraic invariants do not 
diverge faster than their zeroth order piece as the singularity is approached, a natural requirement if one wants the 
first order formalism to provide approximate solution of Einstein's equations that can be consistently interpreted as 
arbitrarily small perturbations of the unperturbed metric. This choice also selects perturbations with finite energy, 
using the energy notion obtained by going to second order perturbation theory [l], Q ■ The singularity of V z at r = r s 
is called a "kinematic" in 0], because it is due to the fact that, as defined, the Zerilli function has a simple pole at 
this point forgeneric smooth gravitational perturbations (see @ and Lemma 1, equation (JTSJ) below.) In the Zerilli 
formulation , the initial value problem (IVP) for linearized gravity around the M < Schwarzschild spacetime 
can then be posed as follows: given ^ z (t = 0,x), ^ z (t — 0, x) defined for x > 0, both satisfying (JTSj) and vanishing as 
x 1 / 2 -or faster- for x — > + , find the unique ^(t,x) obeying the singular equation equation (fi?l)- (fT2")) in the half space 
x > 0, and giving this data for t = 0. 

The purpose of this paper is to solve this rather bizarre IVP. The exterior black hole (r > 2M > 0) Zerilli equation 
is entirely free of difficulties, it is a 1 + 1 wave equation in a complete Minkowskian space (the horizon lying at the 
tortoise coordinate value x — — oo), with a smooth potential, and initial data can be evolved by Ti mode expansion. 
The difficulties for the M < case cannot be overcome by introducing alternative radial variables or integrating 
factors, which can be easily seen to merely move the singularity from the coefficients of the differential equation to 
the measure that makes its radial piece self adjoint. Solving the IVP for M < allows us to complete the proof in 
that the negative mass Schwarzschild spacetime is unstable under linear gravitational perturbations, as part of 
a program to study the linear stability of the most notable nakedly singular solutions of Einstein's equation [5|, i|, 
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in connection to cosmic censorship. Unstable (exponentially growing in time) modes were not only found for the 
negative mass Schwarzschild spacetime 0, but also for the \Q\ > M Reissner Nordstrom and the |J| > M 2 Kerr 
naked singularities [B, @j- The instability for the negative mass Schwarzschild -(A)dS and the negative mass Reissner- 
Nordstrom spacetimes were proved in 0j- The unstable smooth solutions of the M < Schwarzschild linearized 
Einstein equations in [2J, satisfy the desired boundary condition at r = 0, and decay exponentially for large r. It 
is argued in 0] that they can be excited by initial data compactly supported away from r — 0, but this can not be 
proved if we do not know how to evolve initial data. In this paper we show how the IVP for even perturbations 
on a negative mass Schwarzschild spacetime can be solved by using the technique of intertwining operators (see 
Q and references therein). An intertwining operator X = d/dx + g(x) is constructed such that for regular metric 
perturbations ^ := T^> z is smooth and belongs to L 2 ((0, oo), dx). iff satisfies a Zerilli like equation = [d 2 /dt 2 + 7i]W, 
H := —d 2 /dx 2 + V, with a potential V that is free of singularities and such that Ti has a unique self adjoint extension 
in a domain T> C £ 2 ((0, oo), dx), that corresponds precisely to our physically motivated choice of boundary condition 
at x = 0. 

These two key differences with the standard Zerilli approach allow us to give a comprehensive answer to the linear 
stability problem of M < Schwarzschild spacetime, as we can show that physically sensible initial data supported 
away from the singularity generically excites the unstable modes found in [2| ■ As is shown in Section IIII1 this is 
not related to the x = boundary conditions; if a perturbation is initially supported away from the singularity, the 
unstable modes are excited before the excitation reaches the singularity. 

In Section[TT]we give a brief account of Zerilli's approach to (even type) gravitational perturbations of Schwarzschild 
spacetime, stressing the problems that arise in the negative mass case. We exhibit the unstable modes found in 0], 
and introduce the new field which is smooth for smooth metric perturbations, no matter the sign of M, and obeys 
an equation free of singularities for any M. The main results of this paper are listed in a theorem, proved in Section 
IIV1 from which the negative mass Schwarzschild spacetime linear instability follows as a corollary. In Section IIIII 
we illustrate, by means of numerical integrations of the linearized equations, how the unstable linear mode found in 
for the negative mass Schwarzschild spacetime is excited by perturbations with different initial data. Section |V| 
summarizes our results. 

II. SCALAR GRAVITATIONAL PERTURBATIONS OF THE SCHWARZSCHILD SPACETIME 

In the Regge- Wheeler gauge [loj , the scalar perturbations for the angular mode (I, m) are described by four functions 
Ho(r,t), Hi(r,t), H2(r,t) and K(r,t), in terms of which the perturbed metric takes the form, 

/ 2A1\ ( 2M\ 1 

ds 2 = - f 1 —J (1 — eH Y t , m ) dt 2 + 2 eHiY e , m dtdr + ( 1 - — j (1 + eH 2 Y Lm ) dr 2 

+r 2 (1 + eKY l>m ) (d8 2 + sm 2 (9)d^ 2 ) (1) 

where Yi >m = Yi^ m (0,<f) are standard spherical harmonics on the sphere. The linearized Einstein equations for the 
metric (TTJ), obtained by disregarding terms of order e 2 or higher, imply Ho(r,t) = H2(r,t), and a set of coupled 
differential equations for H\, if 2 and K. 

Of particular interest to us is the following unstable solution found in [2| for the negative mass case: 



k{t - r) 
2\M\ 



KM _ "oy 

A(A + 1)[2(A + l)r - 6M] r(r - 2M) t ' 1 \ k(t-r) 

H\(rA) = —H2\t.r) = ~ exp — ; ; — 

V ' ' y,) 36M 2 P [ 2\M\ 



(2) 



where 



and 



(I -1)1(1 + !)(* + 2) 
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(£-!)(£ + 2) 



(4) 



The above solution has the following properties (see Section 7 in 0]): (i) it is exponentially growing in time, (ii) it 
is smooth for r > 0, exponentially decaying for large r; (iii) it has a fast decay as r — > + that guarantees that the 
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first order algebraic and differential invariants of the Riemann tensor do not diverge faster than their zeroth order 
piece, a condition of self consistence of the perturbation procedure, (iii) it has finite gravitational energy Eq 

E g = / G%rf<*dE (3) , V = (1 - 2M/r)- 1 / 2 d/dt, £ = d/dt, (5) 

(2) 

where G ab is the second order correction to the Einstein tensor, and S( 3 ) the spacelike hypersurface orthogonal to rf 
(for details see QBEI.) 

As shown below, the evolution of generic initial data with compact support away from the singularity will excite 
these singular modes, which implies that the negative mass Schwarzschild spacetime is linearly unstable. We first 
recall Zerilli's approach to the linearized problem, in order to exhibit the difficulties in dealing with the evolution of 
initial data for the Schwarzschild spacetime in the negative mass case, and develop an alternative approach to the 
linearized problem that allows us to overcome these problems. 



A. Solution of the even mode linearized Einstein equations: Zerilli's approach 



The linearized Einstein's equation for (fTJ) give a coupled system of partial differential equations involving K, Hi 
and Hi [Io| . This system can be decoupled by introducing the Zerilli function # z (t,r), by the replacements Q, 



where A is defined in (J4j) and 



K 


= q(rY 


Hi 


= h(r) 




d ' 


Hi 


dr 




q(r) = 



2M\ d^ z 



dr 



d^ z 



d 2 ^ z 



dt dtdr 



r I 



h(r)* z + r— 



K 



h(r) = 



A(A + l)r 2 + 3XMr + 6M 2 

r 2 (\r + 3M) 
Ar 2 - 3ArM - 3M 2 



(r-2M)(\r + 3M) 
Note that the relations ((6]) can be inverted and give 



*z(r, t) 



r(r - 2M) 



Ho-r 



dK 
dr 



A + l 



(A + l)(Ar + 3M) 

The full set of linearized Einstein's equations then reduce to Zerilli's wave equation 

a 2 * 



K. 



where, 



dt 2 



H 



+ HV Z = 0, 



dx 2 



V 



looks like a quantum Hamiltonian operator with potential 

2M\ X 2 r 2 [(A + l)r + 3M] + 9Af 2 (Ar + M) 



V = 2 1 



r 3 (Ar + 3M) 2 



and x is the "tortoise" coordinate, related to r by 

dx 
dr 



r J 



We will choose the integration constant such that x — at r = 0, then 

r-2M 



2M In 



2M 



(G) 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



(13) 
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FIG. 1: The left panel shows the I = 2 Zerilli potential for M = 3 as a function of x. 
x — — oo. Note that the potential is smooth and positive. The right panel shows the £ - 
naked singularity is located at x = 0. The "kinematic" double pole is at x ~ 0.761. 



The black hole horizon is located at 
: 2 Zerilli potential for M = —2. The 



1. Case M > 0, stability of the Schwarzschild black hole exterior metric 

For M > 0, the exterior static region r > 2M of the Schwarzschild black hole gets mapped under (fT5|) onto 
— oo < x < oo, with the black hole horizon sitting at x = — oo. The potential V in Zerilli's equation is positive definite 
and behaves as V ~ exp(x/(2M)) as x — ► — oo, ~ x~ 2 as x — > oo (see Figure 1). Equation (JSJ indicates that a 
smooth metric perturbation with compact support in the exterior region corresponds to a smooth Zerilli function in 
L 2 (R, dx). The fact that ty° z := ^ z |(t=o, x ) an d ^1 : ~ d/dt^> z \(t=o,x) can be freely chosen, together with ©, takes 
proper account of the constraints among the initial data for H\,H2 and K. To solve the Zerilli wave equation ^ 
from a given initial data (^°, ty®) 6 i 2 (R, dx) ® L 2 (R, dx), we can use that 7i is a self adjoint operator in L 2 (R, dx) 
to expand ^ z and ^ using a complete set of eigenfunctions of H (TiipE = EifiE)- Equation ([9]) then reduces to 
the following ordinary differential equations for a,E(t) '■= J iPe(x)*^ z(t, x)dx: 

a E = -Ea E 

d E (0) =&%:=[ tl>%9° dx (14) 



whose solution is 



a E (0) = a° E := J dx. 



'o^ cos(VEt) + a% E- 1 ' 2 wn(>/Et) , E > 

ais(*) = { + * 6b ,S = (15) 

a° E cosh(V^Ei) + d E (-£0~ 1/2 sinh(-/ = Et) , E < 

Since the Zerilli Hamiltonian 7i is positive definite, we can use the above equations to obtain and L 2 bound for \t at 
time t in terms of its data ' 131 as 



*| 2 dx < 2 (j \^°\ 2 dx + J *° H- l V dxJ , (16) 

where the inverse of H is defined using its spectral decomposition. The detailed analysis in [ijj gives also the following 
uniform bound for the Zerilli function in terms of the initial data 

\^z{t 7 x)\ 2 < J (j*°| 2 + ^W^° + i|^°| 2 +i>° H' 1 * ) dx (17) 



This proves that the exterior, static region of a Schwarzschild black hole is stable. 
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2. Case M < 0, stability of the Schwarzs child naked singularity 



For M < the range of interest is r > (then x > 0), and a number of difficulties arise due to the fact that q and 
h in |(7J| are singular at r = r s := — 3M/A > 0, and that this point belongs to the domain of interest. This implies 
that ty z is singular at r s , as is also evident from equation (jSJ). The kind of singularity in vf^ is characterized in the 
following lemma: 

Lemma 1: If M < 0, a metric perturbation is smooth if and only if its Zerilli function at any fixed time is C°° in 
open sets not containing r s , and admits a Laurent expansion 



= £ c j(T- r ') j > 



Co 



A 2 C_! 



j>-i 



3M(3 + 2A) 



(18) 



If initial data iff® and \E'° is given, both functions satisfying (|18p. the evolution equation will preserve (118[) . i.e., this 
condition will hold at all times. 



Proof: 

containing r s . A straightforward calculation shows that if ip satisfies (|18l) then so does Hip 
this condition will hold at later times if it is satisfied by the initial data □ 



and 



From ([6]) and 0, the metric perturbation will be smooth if and only if h(r)^f z -r , 
(1 — 2tf/r)^ are smooth. Both conditions lead to (TT8")) . added to smoothness in open sets not 



This guarantees that 



In particular, Zerilli functions for smooth metric perturbations are generically not square integrable (no matter 
which measure we use, either dx = dr/(l — 2M/r) or dr) due to the pole in (fig)) . As an example, the smooth M < 
metric perturbation (0 has the singular Zerilli function 



r (r - 2M) A 
2Ar + 6M 



■ exp 



k(t - r) 
2\M\ 



exp 



kt 



2\M\ 



(19) 



Given that *ff z is a singular function of the metric perturbation, it is not a surprise that the coefficients of the 
differential equation it obeys are singular. This explains the second order pole of the Zerilli potential at r s (and the 
name "kinematic" given in [2j to this singularity.) Note that the approach for solving Zerilli's equation in the M > 
case completely breaks down when M < since (i) <iP z £ L 2 ((0,oo),dx) and (ii) V has the kinematic singularity. In 
particular, the associated quantum mechanical problem with Hamiltonian TL and domain x 6 (0, oo) is not relevant 
in this case because of (i), as discussed in detail in Section 7 of 0. Furthermore, since (0 is a wave equation in the 
half space x > 0, we need to specify boundary conditions at x — 0, besides the initial values of ^ and 4'°, to have a 
unique solution. The fact that the potential has a singularity at the boundary, 



V 



-l/(4x 2 



for x — > 



(20) 



implies that there is an infinite number of (formally, i.e., ignoring the kinematic singularity) self-adjoint extensions of 
H := —d 2 /dx 2 + V, parameterized by 8 G S , obtained by demanding that the Zerilli function behaves as 



cos(<9) 



X 

W\ 



1/2 



sin(0) 



x 

\M\ 



1/2 



In 



x 

\M\ 



(21) 



for x > (the terms in square brackets are the leading terms of two linearly independent local solutions of the 
eigenvalue equation = E^> , E shows up at higher orders). Note that both linearly independent local solutions in 
(|2Tj) are square integrable near x = 0, the potential belongs to the "limit circle class" at x = [3]. This issue was 
analyzed in detail in [l[ (see also [l3| and [14j]) where it was concluded that 9 = is a physically motivated choice, 
since it corresponds to finite energy perturbations with first order contributions to the Kretschmann invariant not 
diverging faster than its zeroth order piece. These results were confirmed in (2j, where it was further shown that 
every algebraic and some of the differential invariants made out of the Ricmann tensor share this property with the 
Kretschmann invariant. Given that this guarantees the self consistency of the linearized treatment, we will be restrict 
our attention to the case \I/ Z ~ x 1 ^ 2 from now on. 

The question left open in [|J is how to evolve initial perturbation data in the M < case, since the Ti mode 
expansion technique used for M > does not apply to the M < case. In the next Section we introduce a field iff 
which is smooth for smooth metric perturbation, and evolves according to a wave equation with a smooth potential 
for any sign of M, thus providing a solution to the initial value problem in the negative mass case. 
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B. Solution of the even mode linearized Einstein equations: alternative approach 

As explained above, the quantum mechanical problem associated to 7i is not directly relevant to the gravitational 
perturbation problem when M < 0. Zerilli's function succeeds in reducing the full set of linearized Einstein's 
equations to a single wave equation, however, for M < 0, this function is singular in the relevant r > range. The 
"kinematic" singularity at r s = — 3M/A in JSj) indicates that physically acceptable Zerilli functions have a simple 
pole at r s (Lemma 1). Thus, even if Ti. could be extended to a self adjoint operator in some subspace of L 2 (lSL,dx), 
this space would not be the natural setting for physically acceptable gravitational perturbations, which, because of 
the kinematic singularity, correspond, generically, to functions that are not square integrable. 

In terms of the Zerilli function, the evolution problem in the negative mass case is: given data it",^ both 
satisfying the conditions in Lemma 1, and vanishing as x 1 / 2 when x — + + (i.e., (|2ip with 9 — 0), find \P z for later 
times. The approach of solving this problem by separation of variables in Zerilli's wave equation and expanding by 
Ti modes fails. A satisfactory solution to the evolution problem requires finding a new field ^ that decouples the 
linearized equations obeyed by K,Hi and H 2 -as z does-, and that is a smooth function of K : Hi and H 2 for any 
value of M. In this section we show how this is done. We will state without proof our main result (Theorem below), 
and illustrate for the I = 2 mode, case in which the explicit formulae are relatively simple. We will defer the proof of 
the theorem to section HVl 

Theorem: Let ipo be the solution of Titpo — given in equations (I52|) and (|55|) . Define g := t/V/V'o (a prime denotes 
derivative with respect to x) and the operators 



J 
J 

Let V be the potential in Zerilli's equation. Then: 
(i) 

7 d 2 d 2 



d_ 

dx 
d_ 
dx 



(22) 
(23) 



V 



dt 2 



dx 2 



V 



dt 2 dx 2 

with V = V - 2g' smooth in the relevant domain (r > 2M if M > 0, r > if M < 0). 

(ii) For any value of M, a metric perturbation {T]) with Zerilli function ^> z is smooth if and only if ^ := is 
smooth in the relevant domain. 



(iii) For M < 0, V ~ 3/(4ir 2 ) as x — ► + . As a consequence H 



c):r- 



+ V(x) has a unique self adjoint extension in 



a domain V C L ((0, oo), dx), defined by the boundary condition ^ 



,3/2 



as x -> 0+ (see PJ,[14|]). Moreover, 



for 4" z as in ([2T|l . z £ T> if and only if = 0. Thus 2? is the set of physically relevant perturbation functions 

(iv) Assume that M < and that (\l/°, is an appropriate initial data set, i.e., it satisfies the conditions in Lemma 
1 and the boundary condition 9 — in (f2"Tj) . Note from (iii) that both T^l\ and Tw* z belong to T>. Let if? be 



d 2 



n 



the unique solution in T> for the wave equation 

conditions * | (t=o,^) = 2"*° and d/dti? \( t =o,x)= 
done in equations ((14"]) and (|15p . The Zerilli field at all times is then given by 



if> = on the half space x > 0, subject to the initial 



This solution can be obtained by H mode expansion as is 



1> z (t,x) 




m(t",x)dt" dt' 



ml 



(24) 



Let us clarify some aspects related to the above theorem. Generically, \& z has a pole at r s (Lemma 1), and so does 
g, then the operators 1 and 1 are singular. The singularities cancel in such a way that *f? :— T'f? z is smooth in the 
domain of interest, that is, X removes the singularity in ^5? z . In the same way — 2g subtracts the pole in V to produce 
a smooth V, As an example, for £ = 2 we have [l| 



■00 



r(r 3 + 3Mr 2 - 6A/ 3 ) 
8M 4 (3M + 2r) 



(25) 
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and thus 



V = 



6(r - 2M)(2r 7 + 5Mr 6 - 9M 2 r 5 - 33Af 3 r 4 - 24M 4 r 3 + 36 Ad 5 r 2 + 36A/ 6 r - 36M 7 ) 

r 4 (r 3 + 3Mr 2 -6M 3 ) 2 



(26) 



whose only singular point lies at r ~ 1.2M, which is outside the domain of interest both for positive or negative Ad. 
The £ = 2 unstable mode HI]) for Af < becomes 



if' 



st = e 2(t-r)/|M| r 3 (r-2M) 4 

4r 3 + 12Mr 2 + 24M 3 



S 2t/|M| Tunst 



(27) 



which is also C°° in the domain of interest. In Figure 2 we exhibit V for £ = 2 and Af = 3 (left), and V for £ = 2 and 

M = —2, superposed with the unstable mode ip unsf (right). 

Since the Einstein's linearized equations reduce to the single equation 



d 2 ~ 



* = 0, *eP 



(28) 



and Tt is self adjoint in T>, we can then solve this equation by TL— mode expansion, in the same way as is done with 
^ z when M > 0. This gives an answer to the issue of evolution in the non globally hyperbolic spacetimc in a way 
entirely analogous to that developed in [l5[ , the only difference being that the radial part of the equations dealt with 
in [l5| are positive, essentially self adjoint operators. Note that we can work entirely in "^-space" with no reference 
to the Zerilli function (as will be done in the next section), and that our alternative formulation also works in the 
positive mass case. The usefulness of equation (|24[) lies in the simpler connection that there is between the Zerilli 
field and the perturbed metric elements if i, if 2 and K, equations ([6]). If we want to construct the perturbed metric 
elements from the shortest way seems to be inserting (|24[) in ©. Note that (|24p is not an evolution equation. It 
just tells us how to recover the information lost after applying XX to the Zerilli function (see Section [TV|) . In fact, 
we need to solve first the evolution problem for x), and then use the solution $(t,x) in the integrand in (|2"4"|) to 
obtain the corresponding solution for ^ z (t,x). 

Corollary: The negative mass Schwarzschild solution is unstable. 

Proof: The spectrum of the operator Ji in T> contains the negative eigenvalue — fc 2 (k given in ((31) ) , with eigenvector 

and C is a normalization constant. Let aP_ k2 and oP_ k2 be the 



^unst _ q x^ uns \ where ip unst is given in JT 
projections of *S?o an d ^0 onto this mode, (see equation l[14p). then from (114[) and (fT5)l applied to 



we obtain 



#(t,a:) 



l -k 2 



cosh(fct) 



_ k 2 



sinh(fct) 



^ unst (x) +Q(t,x), 



(29) 



where $(t, x) is a linear combination of modes ipg, with E 7^ —k 2 . Since the above decomposition is orthogonal 

-1 2 



and thus exponentially growing for large t □. 



cosh(fct) 



l -fc 2 



smh(kt) 



Numerical evidence indicates that ifi unst is the only negative eigenvalue of H. If this is the case, then $ is bounded 
and uniformly bounded in a similar way as the positive mass Zerilli function is, equations (|16[) and (jTTJ) . 



III. NUMERICAL INTEGRATION OF THE EVOLUTION EQUATIONS 

Numerical integrations of the wave equation 



dt 2 



dx 2 



V 



* = 



(30) 
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FIG. 2: The left panel shows the potential V , for 1 = 2 and positive mass M — 3, as a function of x. The black hole horizon is 
located at x — —oo. The right panel shows V for I — 2 and M = —2 (solid line), and the unstable mode ipunst given (|27|l (dotted 
line). This mode satisfies Tfapmst = —jurist- 

subject to the boundary condition 'J ~ x 3 / 2 as x — > + were carried out for ^ = 2, M = — 1 using the Maple built 
in integrator for partial differential equations, working in the standard radial coordinate. The boundary condition at 
r = was enforced by imposing Robin type boundary conditions in the form 3^ — d^/dr = at r = 10~ 4 . We also 
set w = at r — 10 and restricted the initial data and evolution time so that this condition is trivially satisfied. In all 
cases we set ^P = for simplicity. We evolved different 1 $>° initial data sets to see how the unstable Q190 mode gets 
excited, and the resulting numerical solution ^ was contrasted with its expected projection a°_ k2 cosh(fct) ip unst (x) 
onto the unstable mode. 

In the case of Figure 3, ^° = exp(— 10(r — 2) 2 )0(r — 0.0002), a step function. If normalized, this function gives 
a projection a°_ fc2 — 0.79 onto the unstable mode. The unstable mode dominates for t > 1.6, and it is noticeable from 




FIG. 3: Left: evolution of initial data centered at r = 2, with a strong projection onto the unstable mode. Right: time evolved 
unstable mode component for this data. 



t = 0. The evolution of two stable wave packets moving oppositely is also evident in the plot. 

The left panel in Figure 4 shows the evolution up to t = 2.7 of the data <F° = exp(-lO(r-4) 2 )0(r-O.OOO2), *° = 0, 
which has a milder projection onto the unstable mode (a 1 ^ — 0.07 when normalized). The right panel contrasts 

\&(t = 0, r), i&(t = 3, r) and the unstable mode properly scaled by the cosh(6) factor. Note that the unstable mode is 
noticeable starting at t ~ 1.5. 

To have a smaller overlap with the unstable mode we use $ (t = 0, r) = — cxp(— 10(r — 4) 2 ) 0(r — 0.0002) sin(10(r — 
4)) (Figure 5, left panel). The graph may mistakenly (see equation (I29|) ) suggest that the unstable mode is not excited 
before the ingoing wave packet reaches the singularity. The right panel in this figure exhibits the evolution of the 
projection of the initial data onto the unstable mode. This mode is initially highly suppressed because a°_ k2 — 6 x 10~ 3 



FIG. 4: Left: evolution of initial data centered at r = 4. Right: $ at t = (thick solid line), and t = 3 (thin solid line), 
contrasted to the evolution of the unstable component at t — 3 (dashed line). 




FIG. 5: Left: evolution of initial data centered at r — 4 and modulated with a sine function. The initial data has a weak 
projection onto the unstable mode. Right: time evolved unstable mode component for this data. 

Finally, the left panel in Figure 6 exhibits the evolution of data which is almost orthogonal (a°_ k2 — —2 x 
10 _10 ,d^ fc 2 = 0) to the unstable mode. The excitation of the unstable mode reaches an amplitude ~ 2 x 10~ 8 
and so is unnoticeable in the displayed time range. The right panel of the figure shows the initial data (dotted line), 
the result of the evolution at t = 1 (thin solid line) and the evolved data at t = 3 (thick solid line). 

IV. PROOF OF THE THEOREM 

The alternative ^> to the Zerilli field as a means of describing the even modes of linear gravitational perturbations 
is suggested by the "intertwining" potential technique in quantum mechanics [8| , whose original motivation is that of 
replacing a one dimensional quantum mechanical Hamiltonian with another Hamiltonian having a more elementary 
potential. We have actually used an intertwiner to replace a potential with a singularity at r s with one free of such 
singularity, with the added benefit that the resulting Hamiltonian has a unique self adjoint extension that happens 
to agree with the boundary condition at r = that is natural to the problem! Intertwiners appear also the context 
of supersymmetric quantum mechanics pairs of supersymmetric Hamiltonians being related by an intertwiner 
constructed using a zero energy wave function. 



10 





FIG. 6: Left: evolution of initial data centered at r = 2, modulated with a sine function and having a negligible projection 
onto the unstable mode. Right: \& at t = (dotted line), t = 1 (thin solid line) and t = 3 (thick solid line). 



A. Intertwining operators 

Consider a two dimensional wave equation with a space dependent potential V 

* = 0, 



d2 d2 t,/ 



and a linear operator X = — g(x) such that 



J 



d 2 



a 2 



(31) 



(32) 



for some potential V"(ar). Since X commutes with d/dt, any solution W of (|3ip gives a -possibly trivial- solution 
$ := 2"\& for the equation 



d 2 



d 2 



dt 2 dx 2 



V(x) 



* = 0. 



(33) 



Separation of variables = exp(iu)t)ip(x) ( & = exp(icut)tp(x)) reduces (l3"H and (|3"3")l to Schrodinger like equations 

m = -2-+V{x) i/> = u 2 i/>, (34) 
" fl 2 „ 1 . 

= -^ + V{x) i> = J 1 i>, (35) 



If we do not specify boundary conditions, there will be two linearly independent solutions of (fM)) for any chosen 

complex oj. Let us denote any two such solutions as ipu\j — 1,2. Note from (f32|) that Xi/'w are (possibly trivial) 
solutions of (|55|l. 

The conditions for the existence of an intertwining operator can be obtained by applying (|32j) to an arbitrary function 
tp, and then isolating terms in tp and tp' := dtp / 'dx (the higher derivative terms cancel out). The coefficient of tp' gives 



V = V - 2g'. 



(36) 



Adding the condition from the tp coefficient gives (g' + g 2 — V)' = 0, i.e. g' + g 2 = V — lo 2 for some constant cu . This 
last condition is more transparent if written in terms of tp Uc :— exp(J g(x')dx'), which satisfies tp'^ /tp u)o = g and 



dx 2 



V 



IpUo = 0J o tp u 



(37) 



From this follows 0, 
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Lemma 2: From any solution of ([37)) it is possible to construct an intertwining operator X = ^ — g(x) by choosing 
g = ipljtp^. This gives V = V - 2g' in p2]). 

Lemma 2 collects the results we need from ||, but we need to elaborate further on these results to get some 

(j=X) (7 — 2) 

information about the possible ways to invert the effect of X. To fix the notation, let ipu a — j an d ipu a be a 
linearly independent solutions of (|37[) . The kernel of X is the span of since o = xi> = y - ^ =1) 7^ =1 V 

implies that ip is proportional to . The form of an intertwiner X = -j^ — h(x) satisfying 



,2 



2 



(38) 



can be guessed from Lemma 2 by noting that, since V — 2h' = V = V + 2g', the only possible way back to V is that 
n(l/ipi j r 1] ) oc (1/^ =1) ). That this is actually the case can be checked by a direct calculation using our previous 
results, from where we obtain n(l/^L{~ 1] ) = uj 2 /^~ 1] . We will set ^~ 2} := l/^L j ~ 1] and choose V& _1) such that 
I^w^ = V'i'l -1 ''- It follows that X := +g{x) satisfies (|3"5)l . and a simple calculation shows that XXtp = (uj 2 — H)ip, 
i.e., the non trivial kernels of I and X combine in such a way that the kernel of XX is the two dimensional ui 2 eigenspace 
of Ti. Note that we have shown that we can label the solutions of (j3"T|) and its hat version as such that 



^: 


=i) 


= Mr 1} 






=2) 


= ^r 2) = 






= 1) 


= Mr 2) = 


= 



0=1) 



(39) 



We have proved the following 



Lemma 3: The kernel of X = #; — ipH l " /ipu ^ is the subspace spanned by ipw a 1 . If X = d/dx + ip^ /ipu ^ 
then (1551) holds, also 



XX=(co 2 -H), (40) 
and the solutions of TL4>u = w 2 4>u and Hip^ = lo 2 ^^ can be labeled such that the equations (|3"9"|) hold. 



In the supersymmetric quantum mechanics context, lo = and ip^ satisfies appropriate boundary conditions 
to make it an eigenfunction of TL. Moreover, it corresponds to the lowest eigenvalue of H. In this case Ti and Ti are 
isospectral, except for w 2 = 0, which is missing in the spectrum of Ti. Equations (|39p then leads to the situation 
depicted in Fig 2.1 in [9(. In the above construction, however, we do not require any specific boundary condition on 

the function ip^^ used to construct the intertwining operator. 

The intertwining operator (|32p will be useful whenever V is simpler than V. However, information is lost when 
solving (|33p instead of (|3ip. and we need to know how to recover it. This problem is addressed in the Lemma below. 

Lemma 4: Assume ^(t, x) satisfies the wave equation (13"TT) with initial conditions ^(0, x) =: ^>°(x) and d^/dt(0, x) =: 
<F°(x). Let * :=XV, *° := and *° :=X^°(x), then: 

(i) ^> satisfies the wave equation (j3"3"|) with initial conditions \l/(0, x) — ^t and d^/dt(0,x) = Vl/ . 

(ii) If lu 0, *i>(t,x) can be obtained from ^>(t,x) by means of 

V(t,x) = cos(Lu o t)^ + i ' m<yUJot ^ <£°+ — (sm(uj t) [ cos(uj t')X^(t',x)dt' - cos(uj t) [ sm( Lj t')X^ ( t' , x)dt' 



(41) 



For oj n = we have 



*(t,x)= I {f XV(t",x)dt")dt' + ti> + #° (42) 
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Proof: (i) is trivial. To prove (ii) note from Lemma 3, equation (|40]) . that 

= XX* = (cj 2 - H)V = {u) 2 + d 2 /dt 2 )^, (43) 

where we have used that ^ satisfies (j3"Tj) in the last equality. The solution of (|4"3"|) , regarded as a differential equation 
in t on 'J', is 

x) = cos(to t)F (x) + sm(uj t)K(x) + — ( sin(u> t) [ cos(uj t')X^(t' , x) dt' - cos(uj t) [ sin(u t')X^(t' , x)dt' 

Wo V JO Jo 

(44) 

if lu 2 ^ 0, and 

#(*, x) = f jf I*(t", x)dt"^j dt' + tR[x) + Q{x) (45) 

if lu = 0. The unknown functions of x, F and K (Q and i?), are "integration constants" of (|4"3"|) , they contain 
the information about '5 that we have lost when applying XX. Fortunately, this information is just the initial 
conditions, since it can readily be seen that F(x) = ^(O,^) = *ff°(x) and u> a K(x) = d*f?/dt(0,x) = ^°(x) 
{Q{x) = V(0,x) = V°(x), R(x) = d$/dt(0,x) = . This gives gTJ from gU), and © from g5]) □ 



B. Intertwining operator for the negative mass Zerilli equation 

Let H. be the Zerilli Hamiltonian, and assume an intcrtwincr is constructed using a solution of Tiip^ — Lo 2 ijj u . Since 
generic solutions of this equation behave as (|18[) (Lemma 2 in [2(), there is a chance that the transformed potential 
P6p be nonsingular at r s , the singularity of V being removed by —2g', and this may well be a consequence of 

$ = K - (46) 

being a smooth function of the perturbed metric. All these expectations turn out to be right, at least, if we use the 
generalization to arbitrary harmonic number I of the solution of 

Hipo = (47) 

found in Q for I = 2. We will first prove the smoothness of V, then that of <]/ 



1. Smoothness of V 

Given V of the form (fTTjl . M < 0, the transformed potential is 

V = V - 2(^M>)'- 



(48) 



Let us first consider the behaviour of V at the kinematic singularity r = r s . Using the fact that ipo is a solution of 
(|47|) . and turning to r (instead of x) derivatives, we find, 



V = 



2(r - 2M) 2 /(hPo 
r 2 ipQ V dr 



-V(r) 



Now, if ip(r) is any solution of 1[4"7|). 



_! , a Q (£ + 2) 2 (£-l) 



V,(r) = a (r r s )^ + + + ^ ~ ^ + ° ^ ~ ^ 



(49) 



(50) 



where ao, and 0,3 are arbitrary constants. Replacing in (|49[) . assuming ao 7^ 0, and expanding in powers of (r — r s ), 
we find, 



V = 



(I 2 + l + 2)(l + 2) 3 (l-l) 3 
216A/ 2 



(^ 2 +£+ !)(£ + 2) 4 (£- l) 4 
648M 3 



+ l+l) 2 q 3 
3a 



(r-r s ) + 0((r-r s ) 2 ) (51) 
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which shows that V is smooth for r — r s , provided ao 7^ (if ao = 0, V has a second order pole at r s .) We consider 
therefore, -0 O of the form, 



■00 



X(r) 



6M + r(£ + 2)(£- 1) 

with x smooth in r > 0. Replacing (|52|) in (|47|) . we find that x satisfies, 

d 2 X [6M 2 + 2rA(3M - r)] d X [6Af 2 + 2rA(3M + Ar)] 
df2 + r(r-2Af)(3A/ + Ar) dr" ~ r 2 (r - 2A/)(3A/ + Ar) 



x = o, 



then 



V 



2(r - 2AT) 2 /d* 
r 2 x 2 V 



4(r - 2M ) 2 A dx (r - 2Af)(6A/ 2 + Ar(2Ar + 4Af)) 



r 2 (3Af + Ar)x dr 



r 4 (3Af + Ar) 



(52) 



(53) 



(54) 



is smooth at r = r s if x is smooth. The only remaining possible singularities for r > would correspond to the zeros 
of x for r > 0, since V is smooth except at r s . It turns out that (|53"|) admits, for every £ > 2, a polynomial solution 
of the form, 



1+2 



(n - 2) [(n - 4)£(£ + 1) + n - 1] T(£ + n - 1) „ 



2 n T {n) 2 T(£ - n + 3) (-M) r 
which, for i = 2 reduces to the solution (|2"5|) found in [l[, 



3r 



3r 3 



2M 4Af 3 4M 4 

which is positive for AT < and r > 0. Similarly, for M < 0, and £ > 3 we have, 



(55) 



(56) 



x( r ) = 



lAfl 



£(£ + 2)(^ 2 -l) r 2 £(^ + 2)(^ 2 -l) r 3 



32 



96 



|Af| s 



^(£ 2 + £ + 4)(£ + 3)(^ 2 -4)(^ 2 - 1) r 4 



6144 



Af 



(57) 



where all the remaining terms, indicated by dots, are non negative for r > 0. The fourth degree polynomial given 
explicitly between the brackets in (|57| is positive for r = and for sufficiently large r. Therefore, it can only have a 
zero if its derivative vanishes at least at one point for r > 0. One can check that for r > there is only one root given 
by, 



TO 



l\M\ (^U{£ 2 -!){£ + 2) -108 - 6^ 



(£-2){£ + 3){£ 2 +£ + A) 

This must correspond to a minimum of the polynomial in r > 0. Replacing r = ro in (I57p we find, 

16(p - 6) ((£ 3 -£)(£ + 2) - 18) ((£ 3 -£)(£ + 2){ P - 18) + : 



xM > 



(^ 2 +^ + 4) 4 (^-2) 4 (£ + 3) 4 



(58) 



(59) 



where p = ^JU(£ 2 -!)(£ + 2)- 108. The right hand side of ([59]) is positive for £ > 3. We conclude that xO) > for 
r > 0. This completes the proof of the smoothness of V. 

The explicit form of V as a function of r for generic I > 2 is very complicated but, fortunately, it is not required for 
the rest of our analysis. In any case, it is possible to obtain several features of V directly from ([54")) . First, since \ is 
a polynomial of degree £ + 2, we find that for large r, 



v ^ + 2) f + 1 W 3 )>o. 



Also, from (|55p . for r — > we have, in general 

F = 12M 2 r- 4 - 2M(£ 2 +£ + 3)r~ 3 + 0(r~ 2 ) = 



£ 2 + £ - 1 



1(1+1) 

4x 2 ' 4|M| 1 / 2 a; 3 / 2 4|M|a; 



+ 0{x- 1 ' 2 ) 



(60) 



(61) 
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Thus the general local solution of the differential equation HipE = Eipg, for x — ► + is of the form, 

i» B = ao (x 3 ' 2 + •••) + b (V 1/2 + ...) (62) 

which is not square integrable near x — unless bo = 0. This last condition can easily be checked to correspond 
precisely to the 8 = boundary condition for the local solutions of oc VP in (|2Tj) . 



2. Smoothness of ip 

We need only check smoothness at r s , and this follows from equations (dHJ) and (fT5|) . which imply that ip = Tip z 
admits a Taylor expansion around r — r s . Of particular relevance is the transformed of Xi\) unst which, from the 
above results, belongs to D and thus is a negative energy eigenfunction of H, which, therefore, has, at least, one 
bound state. This, by the way, implies that V must have a region where it takes negative values, as can be explicitly 
checked for particular values of £, and is illustrated in Figure 2 (right panel). 



C. Intertwining operator for the positive mass Zerilli equation 

1. Smoothness ofV 

The intertwining transformation is equally applicable when M > 0. One can check that equations (|49|) and (|52 |) - ([55|) 
are still valid if M > 0. From (|54[) we see that for r > 2M, the only possible singularities of V correspond to zeros of 
x(r), then we need to prove that x( r ) nas 110 zeros in r > 2M. This can be seen as follows: first we notice that near 
r = 2M, (|53|) has only one regular solution, and this must correspond to the polynomial solution (|55p . Expanding 
this solution in powers of (r — 2M) we find, 

X (r) = ao + (4 f (2 + A 6 + A 3^1 a ° (' " ™) + ^ - 2Mf + O ((r - 2M)>) (63) 

where ao is a constant. This implies that \ and d\/dr are both non vanishing and have the same sign and, therefore, 
X is increasing if x(r = 2M) > 0, decreasing if x( r = 2M) < 0. Then, in order to have a zero for r > 2M, there must 
be a point where dx/dr = 0. But from 1541 we notice that for r > 2M, at any point where dx/dr = we must have Xi 
and d 2 x/dr 2 with the same sign, namely, this corresponds to a minimum for positive x ancl a maximum for negative 
X- But since, e.g., for x(r = 2M) > 0, the function is already increasing, and the condition dx/dr = cannot be 
satisfied for r > 2M, implying that %(r) has no zeros for r > 2Af , and similarly for x( r — 2M) < 0. Thus V(r) is 
regular for r > 2M, it vanishes as (r — 2M) for r — > 2M (see (J54J) ) , and as l/r^ i+2 ^ e+1 ' for large r. In this respect, 
it is similar to the Zerilli potential V(r). One can see, however, that V(r) is not positive definite, (see Figure 2, left 
panel for an example), making the proof of the stability of the exterior region of a Schwarzschild black hole more 
complicated in the context of the formulation. 

2. Smoothness of 
The proof for M < holds also for positive mass. 

D. Proof of the Theorem 

Parts (i), (ii) and (iii) of the Theorem were proved in the two previous subsections. Part (iv) follows from Lemma 
4 and a uniqueness argument: since there is a unique solution in T> of the equation TC^ = with initial condition 
(\j/° ; \j£r°) g J) x T>, and T^ z is such a solution if ^ z solves Zerilli's equation with initial data ('J' , , J ) and boundary 
condition ^ z ~ x 1 / 2 for x ~ 0, it must be ^ = T^ Zl then part (iv) follows from Lemma 4 and the fact that uj = 
for the intertwiner X that we use. 
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V. SUMMARY 



The propagation of gravitational perturbations on a negative mass Schwarzschild background is a subtle problem 
for two reasons. First, this space is not globally hyperbolic. As a consequence, the perturbation equations can be 
reduced to a single 1 + 1 wave equation with a space dependent potential for the so called Zerilli function, restricted 
to a semi-infinite domain x > 0, (t, x) being standard inertial coordinates on two dimensional Minkowski space, x = 
the position of the singularity. This implies that a physically motivated choice of boundary conditions at x = is 
required. There is a unique choice dictated simultaneously by two conditions [l], 0] : (i) that the linearized regime be 
valid in the whole domain, and, in particular, that the invariants made out of the Riemman tensor behave such that 
their first order piece does not diverge faster than their zeroth order piece as the singularity is approached; (ii) that 
the energy of the perturbation, as measured using the second order correction to the Einstein tensor [l[ be finite. 
The second problematic issue with the standard approach is not essential, but related to a choice of variables: the 
Zerilli function <J z is a singular function of the first order metric coefficients. As a consequence, the wave equation 
it obeys has a potential with a "kinematic" singularity, and it is not clear how to evolve initial data, since the usual 
approach of separation of variables leading to a well behaved quantum Hamiltonian operator for the x coordinate 
breaks down. 

We have introduced an alternative diagonalization of the linearized even mode Einstein's equation around a 
Schwarzschild spacetime, using a field which is smooth for regular metric perturbations, regardless the sign of 
the mass M. This field obeys a wave equation with a smooth potential, that can be solved by separation of variables. 
Moreover, the spatial piece of the modified wave equation has a unique self adjoint extension, that naturally selects 
the boundary condition that is physically relevant. 

The connection between the two fields is provided by an intertwining operator, VP — f z — ip' /ipo^z ='■ 2^« 5 similar 
to the operators linking supersymmetric pairs of quantum hamiltonians. We have also shown that, in spite of the fact 
that X has a non trivial kernel, it is possible to evolve the perturbation equations using, at two different steps, the 
initial condition for the Zerilli function. A straightforward application of this formalism allows us to show that the 
unstable mode found in [2[ can actually be excited by initial data compactly supported away from the singularity. 
This closes a gap in our proof in of the linear instability of the negative mass Schwarzschild spacetime. 
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